UNCLASSIFIED 


_ AD  NUMBER _ 

AD870769 

LIMITATION  CHANGES 
TO: 

Approved  for  public  release;  distribution  is 
unlimited. 


FROM: 

Distribution  authorized  to  U.S.  Gov't,  agencies 
and  their  contractors;  Critical  Technology;  JAN 
1970.  Other  requests  shall  be  referred  to  Air 
Force  Technical  Application  Center ,  VELA 
Seismological  Center ,  Washington,  DC  20333. 

This  document  contains  export-controlled 
technical  data . 


_ AUTHORITY 

usaf  ltr,  25  jan  1972 


THIS  PAGE  IS  UNCLASSIFIED 


249 


a 

o 

00 

a 


SURFACE-WAVE  AMPL I TUDE- VERSUS -D I S TANCE 
RELATION  III  THE  WESTERN  UNITED  STATES 


14  January  1970 


Prepared  For 

AIR  FORCE  TECHNICAL  APPLICATIONS  CENTER 
Washington,  D.  C. 


By 


D.  H.  Von  Seggern 
SEISMIC  DATA  LABORATORY 


Under 

Project  VELA  UNIFORM 


Sponsored  By 


D  D  C 

ihe/hz 

JUN  £4  1970 

EtSEtfTO 'J 


ADVANCED  RESEARCH  PROJECTS  AGENCY 
Nuclear  Monitoring  Research  Office 
ARPA  Order  No.  624 


This  document  is  subject  to  special  export  controls  and  each 
transmittal  to  foreign  governments  or  foreign  nationals  may 
be  made  only  with  prior  approval  of  Chief,  AFTAC.  VS c 


SURFACE- WAVE  AMPLITUDE- VERSUS- DISTANCE 
RELATION  IN  THE  WESTERN  UNITED  STATES 


SEISMIC  DATA  LABORATORY 

AFTAC  Project  No.: 

Project  Title: 

ARPA  Order  No . : 

ARPA  Program  Code  No.: 


REPORT  NO.  249 

VELA  T/9706 

Seismic  Data  Laboratory 

624 

9F1 0 


Name  of  Contractor: 


TELEDYNE  INDUSTRIES,  INC. 


Contract  No. : 

Date  of  Contract: 

Amount  of  Contract: 
Contract  Expiration  Date: 
Project  Manager: 


F33657-69-C-0913-PZ01 

2  March  1969 

$  2,000,000 

1  March  1970 

Royal  A.  Hartenberger 
(703)  836-7647 


P.  0.  Box  334,  Alexandria,  Virginia 
AVAILABILITY 

This  document  is  subject  to  special  export  controls  and  each 
transmittal  to  foreign  governments  or  foreign  nationals  may 
be  made  only  with  prior  approval  of  Chief,  AFTAC. 


This  research  was  supported  by  the  Advanced  Research 
Projects  Agency ,  Nuclear  Monitoring  Research  Office,  under 
Project  VELA-UNIFORM  and  accomplished  under  technical  direc¬ 
tion  of  the  Air  Force  Technical  Applications  Center  under 
Contract  F33657-69-C-0913-PZ01. 

Neither  the  Advanced  Research  Projects  Agency  nor  the 
Air  Force  Technical  Applications  Center  will  be  responsible 
information  contained  herein  which  may  have  been  supplied 
bY  other  organizations  or  contractors,  and  this  document  is 
subject  to  later  revision  as  may  be  necessary . 


TABLli  OF  CONTENTS 


Page  No. 


ABSTRACT 

INTRODUCTION  : 
THE  SURFACE-WAVE  MAGNITUDE  FORMULA  3 
CORRECTION  FACTOR  FOR  REGIONAL  DISTANCES  FROM  NTS  5 
Q  FOR  THE  PERIODS  10-16  SECONDS  13 
CONCLUSIONS  16 
REFERENCES  18 
BIBLIOGRAPHY  OF  SHOT  REPORTS  20 


. 


LIST  OF  FIGURES 


Figure  Title  Figure  N 

Geographical  distribution  of  Rayleigh-wave 

amplitude  measurements.  1 

Individual  magnitude  deviations  from  the 

average  for  31  Nevada  Test  Site  explosions.  2 

Gutenberg  and  NTS  "B"  factors.  3 

AUK  -  and  M^g  4 

BOURBON  -  Mg  and  MWTS  5 

BOXCAR  -  Mg  and  M^Tg  6 

BRONZE  -  Mg  and  M^g  7 

CHARCOAL  -  Mg  and  M^Tg  8 

CORDUROY  -  Mg  and  MNTS  9 

CUP  -  Mg  and  M^TS  10 

FORE  -  Mg  and  M^TS  11 

KLICKETAT  -  MQ  and  MNTS  12 

MISSISSIPPI  -  Mg  and  MNTS  13 

PILEDRIVER  -  Mg  and  M^g  14 

SCOTCH  -  Mg  and  MNTg  15 

TAN  -  Mg  and  M^TS  16 

TURF  -  Mg  and  M^TS  17 

Standard  deviations  of  surface-wave  magnitude 
using  NTS  and  Gutenberg  "B"  factors  for  41 
events. 


18 


LIST  OF  FIGURES  (Cont'd.) 


Figure  Title 


Figure  No. 


Standard  deviations  of  body-wave  magnitudes  using 
Evernden’s  "B"  factor  and  of  surface-wave  magni¬ 
tudes  using  NTS  "B"  factor  for  41  events. 

19 

Least-squares  solutions  for  "B"  factor  for 
ESE  profiles  -  95%  confidence  ellipses. 

N  and 

20 

RULISON-  (Rifle,  Colorado)  -  and 

21 

GASBUGGY  (Farmington,  New  Mexico)  -  M„  and 

M 

NTS 

22 

Effect  of  dispersion  only  on  Rayleigh-wave 
maximum  peak-to-peak  amplitude. 

23 

LIST  OF  TABLE 


Table  Title 


Table  No. 


Least-Squares  Solutions  for  Distance- 
Correction  Factors 


1 


ABSTRACT 


For  epicentral  distances  less  than  15°  in  the  Western 
United  States,  a  new  correction  factor  is  proposed  for  use 
in  the  common  Gutenberg  formulation  for  surface-wave  magni¬ 
tude.  The  data  on  which  this  is  based  are  684  Rayleigh- 
wave  amplitudes  from  Nevada  Test  Site  explosions  measured 
visually  on  the  records  of  LRSM  mobile  stations  and  VELA 
observatories.  The  need  for  the  variable  T  (period)  in 
the  magnitude  calculation  is  discounted  on  empirical 
evidence.  Magnitudes  at  distances  less  than  15°  when 
recomputed  using  the  new  correction  factor  are  in  excellent 
agreement  with  teleseismic  magnitudes  and  show  less  scatter 
among  themselves  than  previously.  An  estimate  of  the 
effective  in  the  crust  from  the  data  is  about  130.  Ampli¬ 
tude  losses  should  reflect  other  causes  than  anelasticity , 

and  this  value  is  undoubtedly  much  lower  than  the  real  Q  . 

K 


INTRODUCTION 


The  determination  of  seismic  event  magnitude  has  a  pro¬ 
longed  history  of  debate  and  confusion.  One  part  of  the  prob¬ 
lem  is  to  formulate  or  tabulate  distance-correction  vales  which 
* 

will,  when  applied  to  observed  amplitudes,  result  in  calculated 
magnitudes  that  are  reasonably  invariant  among  stations  recording 
a  particular  event.  In  addition,  there  have  been  various  pro¬ 
posals  concerning  what  and  where  to  actually  measure  on  the  record 
for  magnitude  determination.  The  proliferation  of  magnitude 
formulas  is  documented  in  the  VESIAC  Advisory  Report  (19G4)  and 
further  by  Bath  (1966  and  1969).  The  generation  of  a  large 
number  of  different  formulas  can  mostly  be  explained  by  either 
differences  in  system  response  among  seismographs  around  the 
world  or  by  differences  in  regional  structure  for  propagation 
paths  studied.  The  first  is  of  course  artificial  and  couid  be 
t liminated;  the  second,  however,  could  produce  significantly 

amplitude- versus -distance  relations,  both  for  body  and 
surface  waves,  in  different  regions. 

This  report  is  concerned  only  with  surface-wave  magnitude, 
and  is  intended  to  serve  as  a  complement  to  part  of  the  work  of 
Evernden  (1967)  which  attacked  the  more  difficult  problem  of 
determining  at  regional  distances  consistent  body-wave  magnitudes 
which  were  comparable  to  those  computed  at  teleseismic  locations. 
The  formulation  of  a  distance-correction  factor  has  largely  been 
an  empirical  effort  in  the  past,  and  we  will  continue  this 
approach  while  giving  only  a  limited  discussion  of  the  theoretical 
causes  of  surface-wave  amplitude  diminution  with  distance. 

The  Seismic  Data  Laboratory  has  accumulated  amplitude  readings 
i or  surface  waves  from  over  fifty  unclassified  explosions  at  the 
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Nevada  Test  Site  (hereafter  referred  to  as  NTS)  and  elsewhere. 

The  advantage  of  using  explosion  data  in  this  type  of  study  is 
that  azimuthally-dependent  amplitude  radiation  patterns  of  earth¬ 
quake  source  mechanisms  will  not  contribute  to  scatter  in  the 
data.  There  is,  of  course,  the  possibility  of  tectonic  release 
accompanying  explosions  as  shown  by  Toksoz  et  al.  (1965),  but 
still  the  explosion  data  should  be  more  suitable  for  this  ampli¬ 
tude  study. 
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THE  SURFACE-WAVE  MAGNITUDE  FORMULA 


The  original  Gutenberg  surface-wave  magnitude  formula  is: 

Mg  =  log10  +  1.656  log1C]A  +  1.818  (1) 

where  A  is  the  sum  of  the  maximum  zero-to-peak  horizontal  ampli= 
tudes  in  microns  of  the  Rayleigh  wave  measured  at  a  period  of 
about  twenty  seconds  and  A  is  measured  in  degrees  (Gutenberg, 

1945,  equation  (4)).  Using  maximum  peak-to-peak  amplitude 
instead,  assuming  a  period  of  twenty  seconds,  converting  to 
millimicrons,  changing  to  measured  vertical  amplitude,  Gutenberg's 
formula  becomes: 

Mq  -  log^g  ( Am^ / T )  +  1.66  log^QA  —  0.18  (2) 

which  is  the  formula  stated  by  Geotech  (1964)  and  employed  at 
the  Seismic  Data  Laboratory  continuously  since  1964.  Note 
•that  one  must  set  T  =  20  to  return  to  (1);  however,  in  practice 
T  is  measured  at  the  maximum  recorded  amplitude  and  is  variable 
in  (2).  Hereafter  in  this  report  the  logarithmic  base  of  10  will 
not  be  written  and  is  to  be  understood;  also,  the  my  subscript 
will  be  dropped  and  all  amplitudes  should  be  considered  as 
measured  in  millimicrons.  This  report  will  speak  of  "Gutenberg 
magnitude"  or  "Gutenberg  distance  correction",  and  it  is  to  be 
understood  that  (2)  is  implied  rather  than  (1)  even  though  (2) 
is  not  Gutenberg's  exact  formulation.  The  term  "B  factor"  will 
also  be  used  to  indicate  (b  logA  +  a)  so  that  (2)  can  be  written 
as 


M  =  log  (A/T)  +  B 


(3) 


-3- 


for  the  general  case  where  1.66  logA  -  0.18  is  not  necessarily 
implied  as  the  correction  factor. 

It  becomes  apparent  over  a  period  of  time,  after  individual 
station  magnitudes  calculated  according  to  (2)  for  various  NTS 
explosions  were  plotted  versus  distance,  that  magnitudes  at 
regional  distances  were  low  relative  to  those  at  teleseismic 
distances.  It  was  also  noted  that  the  periods  of  maximum 
motion  were  almost  always  between  10  and  16  seconds  at  less 
than  teleseismic  distances,  and  so  tne  original  requirement  of 
Gutenberg  that  the  measured  amplitude  be  at  a  period  of  twenty 
seconds  was  seldom  satisfied.  The  explanation  for  the  low  mag¬ 
nitudes  was  simply  that  Gutenberg  did  not  use  regional  data  in 
deriving  his  formula  and  did  not  intend  it  to  be  applied  for 
magnitude  determination  at  regional  distances.  Thus,  the  NTS 
data  warranted  formulation  of  a  distance  correction  factor 
applicable  to  regional  distances  in  the  Western  United  States. 

The  Gutenbe:rg  magnitude  as  given  by  (2)  was  assumed  to  be 
reliable  for  teleseismic  distances  since  he  based  it  on  almost 
one  thousand  amplitude  observations  at  distances  greater  than  13°, 
and  a  regional  distance-correction  factor  was  designed  to  force 
regional  magnitudes  to  agree  with  those  at  teleseismic  distance 
for  the  same  events. 


CORRECTION  FACTOR  FOR  REGIONAL  DISTANCES  FROM  NTS 


As  a  data  base  we  had  available  704  Rayleigh-wave  amplitude 

measurements  from  LRSM  mobile  stations  and  VELA  observatories 

made  visually  at  the  Seismic  Data  Laboratory  and  elsewhere  (see 

Bibliography  of  Shot  Reports)  for  47  NTS  explosions  through 

mid-1969.  Rayleigh  waves  were  recorded  between  1.6° and  40° 

within  continental  North  America.  By  preliminary  analysis, 

twenty  measurements  were  found  to  be  either  highly  anomalous 

or  in  fact  erroneous,  and  these  were  deleted  from  the  data  base. 

Figure  1  shows  the  geographical  distribution  of  the  remaining 

684  Rayleigh-wave  measurements.  To  illustrate  effectively 

the  relatively  low  magnitudes  determined  at  regional  distances 

by  equation  (2),  the  quantity  (M..-M.)  is  plotted  versus  dis- 

^■3  9 

tance  in  Figure  2,  where  FL  is  the  average  magnitude  for  the 

jth  event  as  given  by  the  average  of  the  stations  beyond  15° 

distance  and  FI..  is  the  magnitude  given  by  the  ith  station 

“  J 

for  the  jth  event.  If  there  were  less  than  four  stations  beyond 
15°  for  an  event  no  FL  and  thus  no  (FL^-FL)  were  calculated  for 
that  event.  Eighteen  events  comprising  103  of  the  684  data 
points  were  eliminated  by  this  criterion.  Distances  other  than 
15°  were  used,  but  the  division  at  15°  proved  best.  This  point 
in  distance  roughly  coincides  with  the  nominal  division  into 
"regional"  and  " teleseismic"  distances,  and  mention  of  these 
terms  in  this  report  will  imply  the  division  at  15°  distance. 
Figure  2  provides  support  for  the  assumption  that  Gutenberg’s 
formula  gives  consistent  magnitude  values  independent  of  dis¬ 
tance  in  the  teleseismic  range. 

Conclusive  statistical  proof  of  the  inapplicability  of 
Gutenberg’s  distance  correction  factor  at  less  than  15°  can  be 
given  by  taking  the  684  amplitudes  and  fitting  by  least  squares 
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constants  a  and  b  in  the  equation: 


log  (A/T)  =  -  (a  +  b  logA ) 

for  the  two  distance  ranges  0°  -  15°  and  15°  -  40°  and  again 
for  the  entire  distance  range.  A  method  of  testing  whether 
ooth  groups  of  data  could  possibly  be  fitted  by  the  same 
straight  line  is  given  by  Acton  (1959,  p.  81-82)  and  requires 
that  only  the  standard  deviations  from  the  three  lines  be 
calculated.  An  F  statistic  was  calculated  from  which  it  was 
concluded  with  greater  than  99.9%  confidence  that  the  regional 
and  teleseismic  NTS  data  could  not  be  fitted  by  each  other’s 
least-squares  line.  This  implies  that  the  amplitude-distance 
relation  of  regional  observations  is  definitely  different  than 
that  of  the  teleseismic  observations. 

In  determining  a  regional  distance-correction  factor,  the 
effect  of  the  variable  T  in  (2)  should  be  examined;  it  is 
always  less  than  twenty  seconds  for  the  maximum  amplitude  at 
these  distances.  The  precise  period  of  the  cycle  of  maximum 
amplitude  is  usually  difficult  to  assign  at  regional  distances 
on  continental  structures  since  the  wavetrain  is  dispersed  very 
little;  this  is  a  source  of  error  especially  with  LRSM  and  VELA 
data  since  the  relative  system  response  is  rapidly  decreasing 
at  periods  less  than  twenty  seconds.  Therefore,  in  this  report 
two  other  amplitude  measurements  are  considered,  and  distance- 
correction  factors  are  determined  appropriate  to  them.  The 
first  is  simply  A,  the  maximum  recorded  amplitude  itself,  divided 
by  the  system  magnification  at  the  period  corresponding  to  the 
measured  A  value.  The  second  will  be  designated  A  and  is  the 
maximum  recorded  amplitude  divided  by  the  system  magnification 
at  25  seconds  period,  regardless  of  the  actual  observed  period. 


A 

This  A,  then,  is  the  quantity  least  representative  of  true  ground 
motion  of  the  three  proposed  measures.  However,  even  A/T  and  A 
are  not  necessarily  maximum  ground  motions  since  the  system 
response  distorts  the  true  spectral  amplitudes  of  the  signal. 

The  amplitude  measure  for  magnitude  determination  should  not  be 
judged  on  its  physical  basis  but  rather  on  the  consistency  of 
results  obtained  from  it. 

In  determining  a  distance  correction  factor  it  would  be 
advantageous  to  use  all  the  data  points  available.  This  requires 
a  normalization  of  the  observed  amplitudes,  which  represent 
explosions  having  a  magnitude  range  of  over  two  units,  to  reduce 
scatter  of  this  dependent  variable.  Let  us  assume  the  form  of 
the  regional  distance-correction  factor  to  be  the  same  as  in 
Gutenberg’s  magnitude  formula: 

M  =  log  (A/T)  +  b  logA  +  a  (4) 

where  a  and  b  are  constants  to  be  determined.  To  determine 
the  normalized  value  of  amplitude,  we  write  (4)  for  the  jth 
event  and  ith  station  as 

M..  =  log  (A/T)..  +  b  logA..  +  a  (5) 

-LJ  J-J 

and  define 

N  N 

M.  =  1/N  J  H. .  =  1/N  J  [log(A/T)..  +  b  logA..  +  a]  (6) 

J  i=l  J  i=l  lj 

using  as  the  N  stations  only  those  at  teleseismic  distances 
where  Gutenberg’s  magnitude  formula  is  acceptable.  If  for 
some  event  there  are  less  than  four  teleseismic  stations,  no 
Mj  is  computed  and  any  regional  observations  for  that  event  will 


not  be  normalized  and  will  be  excluded  from  the  least-squares 
determination  of  the  regional  distance-correction  factor.  This 
was  the  case  for  18  of  the  original  47  events  at  NTS,  leaving 
330  points  at  A<15°.  So  is  assumed  to  be  the  correct  event 
magnitude,  and  values  of  the  logarithm  of  the  amplitude  will  be 
normalized  to  it.  Using  M.  rather  than  . ,  (5)  can  be  rewritten 
to  effect  this  normalization: 

log  (A/T)..  -  M.  =  -  (a  +  b  logA . . )  (7) 

i]  3  &  lj  w  ' 

or 

M. 

log  [(A/T)ij/10  3]  =  -  (a  +  b  logA^  )  (8) 

This  can  be  represented  by  the  simple  linear  equation 

yk  =  -  (a  +  bxR)  (9) 

where  yk  is  the  logarithm  of  the  normalized  value  of  ampli¬ 
tude  for  observations  at  regional  distances  for  all  stations 
and  all  events.  The  least-squares  solution  for  a  and  b  will 
give  the  most  consistent  magnitude  values  at  regional  distances 
for  an  event  and  will  simultaneously  force  the  regional  magni¬ 
tudes  to  agree  closely  with  teleseismic  magnitudes  determined  by 
Gutenberg's  formula  for  the  same  event.  Another  approach  could 
have  been  to  tie  the  regional  distance  correction  factor  to  the 
Gutenberg  one  so  that 

a  +  bxk  =  1.56  logA  -  0.18 
at  15° ,  or 
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a  +  b  log  15° 


1.77. 


tven  though  this  would  provide  the  pleasing  result  of  making 
the  distance-correction  curve  continuous,  it  is  more  important 
to  seek  the  best  least-squares  solution  cf  a  and  b  and  accept 
some  discontinuity  at  15°.  Since  Gutenberg’s  data  extended 
from  15°  to  130°,  the  correction  factor  determined  by  least- 
squares  for  this  long  range  is  not  necessarily  accurate  at  the 
15°  end  point. 

Least  squares  solutions  of  (9)  were  computed  for  the  330 

points  at  A<15°  using  the  three  amplitude  measures  previously 

described  so  that  y  represented  log  C(A/T)../10  j],  log 
M.  K  ^  M.  i] 

EA^j/10  -1  J  j  and  log  [A^./IO  j].  The  intercept  a  and  slope  b, 

the  9  5-6  confidence  intervals  on  the  slope  and  intercept  assuming 

the  other  quantity  as  known,  and  the  area  of  the  95%  confidence 

ellipse  on  slope  and  intercept  together  are  given  in  Table  1. 

However,  the  difference  in  area  of  the  three  confidence  regions 

is  not  much,  and  the  small  degrees  of  improvement  gained  by 

As 

using  A  or  A  are  not  really  significant.  The  results  show  that 

A 

A,  a  quite  simple  measure,  is  as  good  as  either  A/T  or  A,  which 
have  more  physical  meaning.  Use  of  A  requires  that  the  analyst 
only  pick  the  largest  amplitude  on  the  record  and  divide  by  the 
magnification  given  at  25  seconds  period,  thus  eliminating  possible 
errors  in  reduction  due  to  assigning  a  wrong  period  and  due  to 
additional  mathematical  operations  when  period  is  involved. 

We  will  illustrate  the  improvement  in  regional  surface-wave 
magnitude  determinations  when  the  revised  "B”  factor  for  regional 
distances  is  used  in  (3)  rather  than  the  Gutenberg  "B"  factor  by 
plotting  magnitude  versus  distance  for  several  NTS  explosions. 
First,  the  regional  (hereafter  called  NTS)  ”B"  factor  and  Gutenberg 
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"B"  factor  are  plotted  in  Figure  3.  Note  that  magnitudes  deter¬ 
mined  at  A=15°  with  the  two  different  MB"  factors  would  differ 
by  0.33  magnitude  units.  This  we  regard  as  a  tolerable  discon¬ 
tinuity  in  the  "B"  factor.  Figures  4  through  17  show  typically 
the  improvement  when  the  NTS  "B"  factor  is  applied  at  A<15° 
rather  than  the  Gutenberg  "B"  factor.  Not  only  is  scatter 
somewhat  reduced,  but  regional  magnitudes  become  compatible  with 
'teleseismic  magnitudes.  This  latter  result  is  very  important 
when  magnitudes  of  small  and  large  events  are  determined  because 
the  smaller  events  are  naturally  recorded  at  only  shorter  distances, 
and  it  is  desirable  that  these  regional  recordings  give  a  magni¬ 
tude  value  equal  to  that  which  would  result  had  stations  at  tele- 
seismic  distances  been  able  to  record  the  event.  A  histogram 
in  Figure  18  has  been  prepared  to  show  the  reduction  in  standard 
deviations  of  surface-wave  magnitude  as  determined  from  all 
recordings  for  those  events  having  at  least  two  LR  observations. 
Evidently  the  use  of  the  NTS  revised  "B"  factor  at  less  than  15° 
produces  a  better  magnitude  estimate  for  an  event  than  the  use  of 
the  Gutenberg  MBM  factor  at  these  distances.  Also,  in  Figure  19 
the  markedly  better  consistency  of  surface-wave  magnitudes  for 
the  41  events  compared  to  the  body-wave  magnitudes  for  these  same 
events  using  Evernden’s  (1967)  corrections  is  evident.  Surface- 
waves  at  regional  distances  are  then  definitely  more  predictable 
than  body  waves. 

We  can  investiage  whether  the  "B"  factor  for  NTS  events  is 
dependent  upon  travel  path  by  restricting  amplitude  observations 
to  certain  azimuths.  The  number  of  observations  in  the  two 
sectors  of  340°--20°  and  110°-130°  epicenter-station  azimuth  are 
sufficient  to  provide  reasonable  confidence  limits  on  a  and  b 
when  solving  (8)  by  least  squares  (using  A  rather  than  A/T).  We 
term  these  two  sectors  the  "N"  and  the  "ESE"  profiles,  and  along 
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them  there  were  57  and  74  observations,  respectively.  Figure 
20  shows  the  results  at  the  95%  confidence  level  where  the 
joint  confidence  ellipses  on  a  and  b  overlap  somewhat.  Also 
shown  is  the  95%  confidence  ellipse  for  a  and  b  using  all  the 
330  observations.  We  can  state  with  a  fairly  high  degree  of 
confidence  that  the  rate  of  diminution  of  amplitude  along  the 
N  profile  is  greater  than  along  the  ESE  profile.  Furthermore, 
there  is  reason  to  believe,  although  with  somewhat  less  confi¬ 
dence,  that  diminution  of  amplitude  throughout  the  rest  of  the 
Western  United  States  is  greater  than  along  either  of  these 
profiles  since  the  result  for  all  observations  is  weighted 
heavily  (60%)  by  observations  outside  these  profiles.  These 
differences  in  diminution  rates  can  be  attributed  to  regional 
tectonic  nature,  and  this  will  be  discussed  later.  However,  "B" 
factors  determined  for  the  two  profiles  and  for  the  entire 
group  are  not  more  than  .11  magnitude  units  apart  at  any  dis¬ 
tance  from  1°  to  15°,  and  so  for  practical  application  in  the 
Western  United  States,  the  NTS  "BM  factor  determined  from  the 
entire  group  of  data  would  be  sufficiently  accurate  even  along 
the  N  and  ESE  profile. 

To  test  the  applicability  of  the  NTS  MBM  factor  for  events 
located  elsewhere  in  the  Western  United  States,  Figures  21  and 
22  were  prepared  in  the  manner  of  Figures  4  through  17  for  two 
recent  shots,  RULISON  and  GASBUGGY,  respectively.  RULISON  was 
detonated  near  Rifle,  Colorado^  and  the  data  was  taken  from 
preliminary  analysis  at  the  Seismic  Data  Laboratory.  Magnitude 
data  for  the  GASBUGGY  event,  near  Farmington,  New  Mexico,  was 
taken  from  Rasmussen  and  Lande  (1968).  Figure  21  shows  that 
definite  improvement  is  made  by  use  of  the  NTS  ”BM  factor  for 
the  RULISON  event  even  though  it  is  about  700  kilometers  to  the 
east  of  the  Nevada  Test  Site  and  in  a  different  tectonic  region. 


The  GASBUGGY  event,  Figure  22,  reveals  some  improvement 
with  the  use  of  the  NTS  "B"  factor,  but  regional  magnitudes 
are  now  overestimated  relative  to  teleseismic  ones;  however, 
the  teleseismic  magnitudes  for  this  event  are  few  and  are 
very  scattered  so  that  the  average  teleseismic  magnitude  has 
quite  wide  90%  confidence  limits  of  +_  0.25  magnitude  units. 
There  were  no  Unite d  States  exp lo s ions  outside  of  the  Wes  tern 
United  States  on  which  to  test  the  NTS  "B"  factor  mainly 
because  of  insufficient  teleseismic  recordings. 


-12- 


Q  FOR  THE  PERIODS  10-16  SECONDS 


The  diminution  of  surface-wave  amplitude  with  increasing 
epicentral  distance  can  be  wholly  attributed  to  these  causes: 

1)  geometrical  spreading,  2)  defocusing,  3)  dispersion,  4)  ab¬ 
sorption,  and  5)  various  reflection,  scattering,  energy 
conversion,  and  acoustic  amplification  processes  inherent  wave 
propagation  through  an  inhomogeneous  medium  with  non-parallel 
layering.  Of  these  causes,  only  the  first,  geometrical  spread¬ 
ing,  is  uniform  and  invariant  over  the  earth.  Harkrider  (1964, 
equation  (85))  shows  that  the  distance-dependence  of  the  verti¬ 
cal  Rayleigh-wave  displacement  from  an  explosive  source  at  the 
^ flee  of  a  multilayered  media  is  expressed  by  the  zero— order 
Hankel  function  of  the  second  kind  [H^2)  (kRr)]  wh^re  kR  is  the 
wavenumber  and  r  the  epicentral  distance.  Whenever  kRr  is  greater 
^1/ 2ak°U^  unity’  function  can  be  closely  approximated  by 

r  times  a  constant,  and  so  cylindrical-wave  spreading  closely 

approximates  the  Rayleigh-wave  spreading  at  distances  as  close 
as  100  kilometers.  The  sphericity  correction  for  Rayleigh  waves 
out  to  15°  is  negligible.  Thus  considering  geometrical  spread¬ 
ing  only,  amplitude  will  be  approximately  proportional  to  A-1/^2 
for  the  distance  range  of  this  study.  The  second  cause  of  ampli¬ 
tude  diminution,  defocusing  as  discussed  by  McGarr  (1969),  can 
probably  be  disregarded  in  this  study  because  the  extensive  coverage 
over  the  Western  United  States  of  the  data  used  would  tend  to 
average  out  defocusing  effects  with  similar  focusing  effects. 

To  study  the  effect  of  dispersion  on  the  amplitude-distance 
relation,  several  phase-velocity  dispersion  curves  were  selected 
from  the  collection  of  Brune  (1969,  and  extrapolated  if  necessary; 
and  Rayleigh-wave  signals  comprising  the  period  range  5-100  sec¬ 
onds  were  synthesized  in  the  manner  presented  by  Sato  (1960)  at 
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distance  increments  of  200  km  out  to  1600  km  using  as  a  pre¬ 
sumed  source  spectra  the  Fourier  amplitude  spectra  of  the 
CORDUROY  event  recorded  at  KN-UT,  an  epicentral  distance  of 
about  300  kilometers.  A  time  window  of  the  KN-UT  record 
corresponding  to  Rayleigh-wave  group  velocities  of  4.1  and 
2.8  kilometers  per  second  was  used  in  the  Fourier  transfor¬ 
mation;  the  end  points  of  the  time  window  are  not  critical 
since  the  maximum  recorded  amplitude  occurs  at  some  inter¬ 
mediate  group  velocity.  The  results  of  peak-to-peak  maximum 
amplitude  measurements  on  the  synthesized  signals  for  four 
different  crust  and  upper-mantle  structures  is  presented  in 
Figure  23  using  Brune ' s  terminology.  The  log-linear  plot 
shows  that  the  data  follows  a  relation  of  the  form  A=8e~Yr. 

The  value  of  the  exponent  can  be  determined  only  crudely  since 
the  data  used  in  this  study  crosses  several  varying  tectonic 
regions.  We  select  though  a  value  of  the  exponent  which  is 
midway  representative  ok  the  "Basin-range"  and  the  "Mid¬ 
continent"  structures,  And  this  is  about  -0.00018r. 

We  can  now  correct  lour  330  regional  time-domain  ampli¬ 
tudes  for  known  geometrical  spreading  and  dispersion  effects 
by  the  relation 


a-*  _  a  .1/2  .  0  2  0  A  .  . 

A .  .  =  A .  .  A .  .  e  in 

13  13  ID  J 

where  A^\  is  the  corrected  amplitude  corresponding  to  the 
observed  amplitude  A^j  at  station  i  for  event  j.  These 
corrected  values  can  be  normalized  to  average  teleseismic 
event  magnitude  Mj  as  done  previously  and  then  used  in  the 
common  equation  for  attenuation, 


[At./10Mj]  =  A  e"aA 
i]  o 
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to  estimate  by  least  squares  the  absorption  coefficient  a  over 
the  Western  United  States.  The  value  of  a  was  found  to  be  .067 
(degrees)  with  95%  confidence  limits  of  +  .024,  and  since 
almost  all  of  the  amplitudes  used  in  this  fit  were  picked  at 
periods  of  10  to  16  seconds,  this  value  will  apply  only  to  this 
limited  band.  In  more  familiar  units,  a  =  .00060  km-1  and  the 
quality  factor  QR  follows  from  the  relation 

Qr  =  fr/aUT 

where  U  is  the  group  velocity.  For  our  data  QR=sl34  with  95%  con¬ 
fidence  bounds  of  98  and  208.  This  is  to  be  regarded  as  an  "ef  f  ec- 
tlVe"  Qr  Slnce  we  have  not  separated  out  the  true  absorption 
effect  from  other  causes  of  amplitude  diminution  such  as  reflec¬ 
tion,  scattering,  and  mode  conversion  of  Rayleigh  energy  in  the 
Western  United  States.  The  complex  topographic  and  tectonic 
haracter  of  the  area  covered  by  the  data  suggests  that  these 

latter  phenomena  may  signif icatnly  lower  the  effective  Q 

XR 
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CONCLUSIONS 


Using  330  measurements  of  surface-wave  amplitude  at 
regional  distances  along  with  teleseismic  measurements,  a 
"B"  factor  for  the  Gutenberg  magnitude  formula  has  been 
determined  which  is  applicable  to  the  Western  United  States 
At  less  than  15° ,  the  relation  is 

M  =  log(A/T)  +  1.16  logA  +  0.74. 

Values  of  M  determined  at  A<15°  by  using  A  or  A  show  some¬ 
what  less  variance,  and  use  of  these  measures  is  preferable 
in  practice  since  reduction  to  ground  motion  from  the  film 
record  involves  fewer  steps  and  removes  a  variable  from  the 
magnitude  determination. 

The  differences  in  attenuation  along  particular  paths 
around  NTS  is  insignificant  in  regard  to  magnitude  values, 
and  the  data  used  does  not  warrant  determination  of  path- 
dependent  "B"  factors  in  the  Western  Unites  States. 

The  effective  (134)  in  the  Western  United  States  is 
quite  low  for  the  10-16  second  period  range.  Since  the  Rayleigh 
waves  for  these  periods  are  primarily  contained  in  the  crust, 
this  value  is  representative  of  the  crustal  layers,  and  as 
such  is  lower  than  that  from  the  MM8  model  of  Anderson  et  al. 
(1965)  or  the  "high-frequency"  model  of  Tsai  end  Aki  (1969). 

It  must  be  pointed  out  though  that  these  models  were  not  con¬ 
structed  using  data  with  periods  as  short  as  used  in  this 
report.  However,  a  value  of  Qg~450  for  shear  waves  in  the  crust 
was  determined  by  Press  (1964)  using  the  L  phase  from  NTS 
explosions.  Further,  a  value  of  Q^IOOO  for  compressional 
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wavej  in  "the  crust  was  determined  by  Archambeau  et  al.  (1969  ) 
using  the  phase  from  explosions  in  Nevada.  Since  theoretically 
(Ancer >on  et  al. ,  1965)  QR  should  be  about  5%  greater  than 
and  about  one-half  Q^,  our  value  of  QR  =  134  indicates  definitely 
that  diminution  of  Rayleigh-wave  amplitude  in  the  Western  United 
states  is  due  to  causes  other  than  anelasticity  if  we  assume 
that  LR  waves  are  sampling  much  the  same  material  as  these  body- 
wave  phases  and  that  Q  is  independent  of  period  over  the  interval 
of  one  to  sixteen  seconds. 

The  difference  in  attenuation  rates  along  the  and  ESE  pro¬ 
files  from  NTS  suggests  that  the  effective  QR  is  path-dependent 
in  the  Western  United  States.  Whether  this  represents  real 
changes  in  the  anelasticity  with  region  or  other  contributing 
factors  to  amplitude  diminution  is  a  formidable  problem  requiring 
more  sophisticated  analysis  than  undertaken  here. 
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Figure  2.  Individual  magnitude  deviations  from  the  average 
for  31  Nevada  Test  Site  explosions. 
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Figure  18.  Standard  deviations  of  surface-wave  magnitude 
using  NTS  and  Gutenberg  "B"  factors  for  41 
events . 
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igure  19.  Standard  deviations  of  body-wave  magnitudes 
using  Evernden's  "B"  factor  and  of  surface- 
wave  magnitudes  using  NTS  "B"  factor  for  41 
events . 
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tances  less  than  15°~w>ien  recomputed 
ctor  are  in  excellent  agreement  with 
show  less  scatter  among  themselves  than 
the  effective  Qn  in  the  crust  from  the 
de  losses  shoula  reflect  other  causes 
value  is  undoubtedly  much  lower  than 
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